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Abstract. New expansions for some functions related to tlie Zeta function 
in terms of the Pochhammer's polynomials are given (coefficients fe^, d^, 

and dfc). In some formal limit our expansion 6fc obtained via the alternating 
series gives the regularized expansion of Maslanka for the Zeta function. The 
real and the imaginary part of the function on the critical line is obtained with 
a good accuracy up to 3{s) = t < 35. 

Then, we give the expansion (coefficient dk) for the derivative of ln((s — 
1)C(«)). The critical function of the derivative, whose bounded values for 
SH(s) > I at large values of k should ensure the truth of the Riemann Hy- 
pothesis (RH), is obtained either by means of the primes or by means of the 
zeros (trivial and non-trivial) of the Zeta function. In a numerical experiment 
performed up to high values of k i.e. up to fe = 10^^ we obtain a very good 
agreement between the two functions, with the emergence of twelve oscillations 
with stable amplitude. 

For a special case of values of the two parameters entering in the general 
Pochhammer's expansion it is argued that the bound on the critical function 
should be given by the Euler constant gamma. 



1. Introduction 

Lately there has been new interest in the study of the expansion of the Zeta 
function via the Pochhammer's polynomials. This is related to the original idea 
of Riesz [1] and of Hardy-Littlcwood [2] at the beginning of the last century. In a 
pioneering work [3] Maslanka obtained a regularized expansion for the Zeta function 
(with coefficients A^) and Baez-Duarte for the expansion of the reciprocal of the 
Zeta function (with coefficients c^) for the Riesz case [4, 5]. Other cases of interest 
have also recently been studied [6, 7, 8, 9, 10]. As pointed out in [5], the discrete 
version by means of the Pochhammer's polynomials Pfc(s), where s = a + it\s the 
complex variable and k is an integer, has advantages especially in the context of 
numerical experiments in connection with some "kind of verification" in the direction 
to believe that the RH may be true. 

In this work we first derive a new expansion for the Zeta function in terms of the 
Pochhammer's polynomials via the alternating series (with new coefficients bk)- In 
some formal limit, a connection with the expansion of Maslanka is also obtained in 
Section 2. Our expansion is then studied numerically on the critical line where a 
good agreement with the real function is obtained up to Di(s) = i < 35, with the 
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emergence of the first few low zeros. After this value of t, a divergence possibly of 
numerical nature set on. 

In Section 3 we then obtain the expansion for the function ln((l — 2^~*)C(s)) 
(with new coefficients dk) as well for the derivative of ln((s — l)C(s)) (with new 
coefficients dk) in terms of the two parameters a and /3, already introduced in 
our previous works [11, 12, 13] Then the critical function for the derivative (whose 
boundedness at large k would ensure the truth of the RH) is then obtained either 
with the primes or with the trivial and non-trivial zeros of the Zeta function. 

In a numerical experiment for the special case a = | and ft = 4 up to high values 
of k, i.e. k = 10^^, the results for the two functions are in very good agreement, 
both with the emergence of the same twelve oscillations of stable amplitude of about 
0.01 (Section 4). 

Finally, in the limit of large ft and a = 1, it is argued that an upper bound to 
the critical function should be given by the Euler constant gamma (Section 5). 



2. Zeta function representation via the alternating series 

In this section we derive a formula for (1 — 2^~*)C(s) similar to the one of 
Maslanka [3] for (s - l)C(s) and of Baez-Duarte [4, 5] for [C(s)]~^ 

Here the starting series is convergent for fH(s) = a > and the formula is 
obtained still in terms of the so called Pochhammer's polynomials of degree k, in 
the complex variable s = a + it. 

k 

(2.1) Pk{s) = n " -) Vfc e N* and Po{s) = 1 

We will also use a family of functions with two parameters {a and ft) as considered 
already in our recent works [12, 11, 13]. Since the alternating series is given by: 

(2.2) (1 - 2i-«) as) = J2 ^^(«) = ^>^ 

1 ^ 
n— 1 

we have using the trick as in [4] that: 

(i-2-«)C(^) =E^(i-(i-^))'^ 

n=l 

= E^^E(-if(i-^r(T) 

n=l k=0 

Since 

(-i)'(T) =^(^ + i-i)---(^ + i-'^) 

we obtain: 

(i-2i-)C(s) = E^/c(^ + i)E ^^(1-;^)' 
= E^'fe(^ + i)E(-irO E 

fe=0 j=0 n=l 
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Since from (2.2) 



oo , sn— 1 



(-1) 



n=l 



(l_2i-("+/3.))C(a + /3j) = ^ 

n=l 

substitution in (2.3) gives: 

oo k / 7 \ 

(2.4) (l-2i-)C(3) = E^'^(^+l)E(-l)n J (l-2^-^"+''^0^("+/3i) 
With the definition 

k 

(2.5) fefe := (-1)' 0) (l - 2i-("+«)) C(a + 
(2.4) becomes: 



(2.6) (1 - 2^-) C(.) = bkPkC—^ + 1) 

k=o ^ 

where Po(^ + 1) = 1 and 6o = (1 - 2^-")C(a). 

The scries above, is expected to represent (1 — 2^~'')({s) for s in some compact 
subset of the plane as for the Maslanka case [3]. In that case, the central point has 
been investigated and elucidated by Baez-Duarte [14]. Here many choices of a and 
/? are possible. For a = /3 = 2 we have the Riesz case [1] and it is the analogon to 
the regularized version of Maslanka but the representation of the Zeta function is 
not the same. For a = 1 + S {S iO) and /3 = 2 we obtain the Hardy-Littlewood case 
[2] which was also discussed numerically in a different way using other polynomials 
[15]. 

In fact, from Lemma 2.3 of Baez-Duarte [5] which states that at large k: 

(2.7) \Pk{s)\ < Ck-'^^'^ 

where C is a constant depending on |s|, 
we obtain here that: 



We thus suspect and expect that the above series represents (1 — 2^ *)C(s) for 
all 9^(s) > I + (5, (5 > if we assume |6fc| < Dk~'^ with 7 > at large values 

of k and for some constant D. In fact with this assumption we have that: 



(1 - 2i-«) C(s)| < E bkPki^ + 1) < const. E k ^k-i^^+^) 

' " fe=0 



fe=0 



< const. A; ^ + 1^ ' < 00 

k=0 

if m{s) > ^+S. 

For a = P = 2 (Riesz) we should have \bk\ < Dk~i^'^. For the case a = 1 and 
fi = 2 (Hardy-Littlewood) we should have 

\hk\ < Dk-^+\ Another case of interest 
is the one where a = | and 3 = 1. In this case one should have \hk\ < Dkr^^'' . 



Of interest also, is the limiting case of large values of /3, where barely 6^ should 
behave as \b}.\ < D. 
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For a strong argument (a Theorem) in favour of the vaHdity of the Maslanka 
representation of (s— l)^(s) in some regions of the complex plane (compact subsets), 
the reader should consult the works of Baez-Duarte [14] already mentionned and 
it is expected that using the same methods, the proof of (2.6) may be obtained for 
all 9\{s) > \. Here, for our series we limit ourselves to a numerical analysis just 
illustrating the kind of accuracy of some representations. 

Remark. Let us consider the Riesz case a = (3 = 2. We can write: 

(l - e^'-'^'-') as) = ^2^4) E ( ) (l - C(2 + 2j) 

fe=o j=o ^•'^ 

and using the Taylor's expansion of e^, we obtain: 

oo 

(2.8) (s-l)C(s) = ^A,Pfc(|) 

fe=0 

where 

(2.9) A, = i-iy Q (2j + 1) C(2j + 2) 

i.e. the representation obtained originally by a different method by Maslanka in a 
pioneering work [3]. We remark that (2.8) and (2.9) should not be considered as an 
approximation of our formulas (2.5) and (2.6) and vice versa. (2.5), (2.6) and (2.8), 

(2.9) are simply two different representations of functions related to the Riemann 
Zeta function, the first one given by (s — l)C(s)! the second one by (1 — 2^^-'^)({s). 

As an example, for s = a with a in [0,1], both representations give a good 
description of the real function ({a) as may easily be computationally checked. 

We now proceed to obtain a representation of ({s) possibly correct on the critical 
line s = ^ + it, with the help of (2.5) and (2.6), in which we are free to set a = | 
and P = i. Then: 

(2.10) (l - 25-'*) C(2 +it) = Yl ^^P^^'^ + 1) 

fe=o 

where now 

k 

(2.11) h = E (-1)^ Q) (l - 2^-^) C{\ + ij) 

We now check the scries in (2.10) restricting k up to 20 for t < 18 and up to 50 
for t > 18. We compare the result with the exact functions 9^((1 — 2*)(^(s)) and 
- 2*)C(s)), for s = l+it with t up to 40. The plots are given below. We 
obtain a good approximation with the emergence of the first five non-trivial zeros 
located at ti = 14.13472 . . ., ^2 = 21.02204 . . ., ts = 25.01085 . . ., ti = 21.02204 . . ., 
t5 = 32.93505 The numerical results are satisfactory until t = 35. 

This concludes the first part of our work. Below, in the second part we develop 
two new representations of the functions ln((l — 2^~^)C,{s)) and ^ln((s — l)C(s)) 
which may possibly constitute a satisfactory approximation to the exact functions. 
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3. A REPRESENTATION FOR THE LOGARITHM OF THE ZETA FUNCTION AND AN 
ADDITIONAL CRITERION FOR THE TRUTH OF THE RH 

We will start as before but instead of writing C(s) as a sum, i.e. C(s) = 
X] 7^, we will use the Euler product formula to derive a new representation for 

n=l " 

ln((l — 2^^'*)C(s)), which of course should be carefully investigated by means of 
some numerical experiments. Thus: 

(3.1) ln[(l-2i-)C(.)]=ln[(l-2i-) J] Vin(.) > 1 

p prime 
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For any prime p, we have: 



so that introducing the parameters a and f3 as before we have that: 



n=l 



fe=0 

E Pk{^ + i) E ^ E 

fe=0 n=l j=0 

E + 1) E (-!)'( ■) ln(l 

/c=0 j=0 



the same treatment for the function lii(l — 2^ gives: 

oo A; / 7 \ 

in(i - 2^-) = x: PkC-^ + 1) E (-i)u J - 2' 



where are still the Pochhammer's polynomials. 

Finally, the representation of lii((l — 2^~^)<^{s)), we propose is given by: 



\n[{l-2'-'^)as)] = J2dM- 



k=0 



(3 



1) 



(3.2) 
where now: 

(3.3) 4 := E (-1)' ( •) - 2^-("+^^-)) Q{a + m 

Remark. Another formal derivation of the above equations is the following: 



ln[(l - 2^-^) as)] = In 



E 

.n=l 



Supposing now that the right hand side may be given as an unknown series 
E we then have: 



r=l 



E^(i-(i-^))^ 

r=l 



fc=0 r=l 

cso oo 

E^'fe(^ + i)E(-ir(-)E^ 



fc=0 



J=0 



E^fc(^ + i)E(-iy(-)in( E ^ 



fe=o 



(-1)" 



n=l 



which coincide with (3.2) and (3.3), obtained with the Eulcr product formula for 
5H(s) > 1. (3.2) with (3.3), is the new formula possibly representing the logarithm 
of the Zeta function in terms of the two parameters Pochhammer's polynomials. 
To the best of our knowledge the above representation is new and it is our aim to 
carry out some numerical experiments in the sequel in order to support its validity 
also in some compact subset of the critical strip. 
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We now investigate the representation of the derivative of the function ln((s 
l)C(s)): 

Then with C(s) = Y\ i^^-s we obtain: 

p prime 

=-E^ln(l-p-^) = -ET4=.^(l-e-«^"-) 

p p 

p p g=l 

Introducing as above the Pochhammer's polynomials we obtain further: 

CM 



--Ein;^E^(i-(i-^)) 



oo k . oo 

= - E inp E + 1) E (-1)' C') E 

p k=0 i=0 g=l ^ 

oo A' oo / \ 

= E + 1) E (-ir ( ■) E - E Inp 

oo fc ' k d I I ll 

/c=0 i=0 \ g=l p y 

CXI k . [ 

= E + 1) E (-ir 0) ^ - E in(i - ^) 

fc=0 i=0 \ P J 

oo 

= E + E (-i)^C)|^in(n 

oo k 

= E ^^^(^ + 1) E (-irO)^lnC(a + /3j) 

fc=0 j=0 

For -j^, using -^zri= Jo° e~^^*~^''(iA we have similarly: 

= E + 1) E {-iy{';)^^dx 

k=0 j^O 

oo k 

= E ^^^(^ + 1) E (-l)^0re-^("+/3^-i)dA 
i=o 

oo ^ ■ J 

= Em^ + i)E(-i)^C)^^+^ 

fc=o i=o 
= E ^(^ + 1) E (-irO|^ln(a + /?j-l) 

fe=0 j=0 

Thus, along these lines we obtain: 

(3.5) I His - 1) C(s)) = £ 4Pfe(^ + 1) 

fe=o ^ 
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where: 



(3.6) 4 = E (-1)' Q) ^ I'^K" + ^3 - 1) c(« + m 



]=0 

From the formula (7) in [16], where p represents a non-trivial zero of the Zeta 
function, i.e.: 

1 I C (s) _ s_ I V i -u V V J- -u V 1 4_ C (0) 

5-1 "T" C(s) ~ s-1 s-1 p s-p 2n ^ s+2n "'^ C(0) 

/5 p n—1 n=l 

, , OO OO 

~ C(0) -^^Z^p Z^2n^Z^ s-p ^ s+2n 
p n=l p n=l 

Setting C = — 1, this equation applied to s = a + Pj in (3.6) gives: 
4 = E (-l)' 0\^ + /r ( E e-^(«+«-^) + e-^P + E e-^(«+«+2") - e-^^A dX 

3=0 ^ \ \ p n=l J J 

= r E (e-^("-'')(l - + e-\l - ^ts,^ dX+ 

We consider only fc > 0. Now we make the variable change e~'^^ = x and finally 
we obtain: 

4 = ^ lo' (1 - xr'-' E x'^-'dx + i (1 - .)'=+^-^ E x'^-'dx 

p n=l 

= ^ES(^^,fc + i) + i E B{^,k + i) 

p n=l 

where -B(a;, y) = ^^^J^^^'* is the Beta function. Thus for large k we can write: 

(3-7) 4 = ^Er(^)^-^ + ^ Er(-^)fc-^ 

For the critical function [12] corresponding to 9^(s) = a we have an analogous 
expression to the Baez-Duarte formula for the Cfe appearing in the expansion of 
Cis)-' [4, 5]: 

(3.8) ^ " 4 = ^Er(^)^ ' )^ ' =■ Mk) 

On the other hand we can express 4 and then the critical function with a second 
formula: 

(3.9) A = lr(!i^)*--i'- £ ,„,f;-L(i -i;' 

p prime g=l ^ 



a — o- ^ 1 rV 1 l-cr CK-cr x ^ X > 1 / 1 

(3.10) k^d, = -ri—)k^-k^ Y: i-^'Ep^fi-^ 

p prime g=l ^ 



fc 

=: Mk) 
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In fact (see above) the Pochhammer expansion for is: 

fe=0 ^ 

where 



which for large A; behaves as jjT{^j^)k f . Indeed with the substitution e '^^ = x 
we obtain: 

1 1 fc+i-i 1 Q.-1 

Sk = -s X 1^ \l-x)dx=- xi^ \l-x) dx^ -B{——,k+l) 
P Jo P Jo P P 

It is interesting to note that one can express the critical function in terms of the 
zeros of the Zeta function (3.8) or in terms of the primes (3.10). We will investigate 
numerically these two functions for the case a = |, /3 = 4, a ~ ^. 



4. Numerical experiments 

As a test of the goodness of (3.2) we draw in Figure 3 the plots of the function 
ln((l — 2^~'^)C((t)) and of its polynomial representation in the interval a G [—1, 1[. 
Figure 3 shows a good match between them also in the "critical real interval" [0, 1]. 




Figure 3. The function ln((l - 2'^-'^)({a)) [black] and its poly- 
nomial representation [red] 



In the next figures we present the results of the numerical experiment performed 
on our representation (3.5) for the case a = | and P — 4. We calculated the 
critical functions ipi and V'2 for ^{z) = a = \. In our calculations we considered 
only the first 10 non-trivial zeros of the Zeta function, the first 20 trivial ones and 
the first 5'000 primes. Furthermore using the usual substitution x = \ogk, V'l and 
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ip2 become: 



= H Era - 'i)e'^ + Er(i + f + Er(^-+ 

\j=l j=l n=l 



50 



g=l 



5000 

primes 



where tj is the imaginary part of the j-th non-trivial zero. 

We argue ^2 should approach -01. The convergence is surprising. The computa- 
tions presented in Figure 4 and Figure 5 indicate that the qualitative and quanti- 
tative agreement between the two functions is very good in the range 2.5 < a; < 30 
(15 < < 1.068 X 10"). 




log k 



-0.01 



Figure 4. The critical function ipi calculated with the zeros of 
the Zeta function 




log k 



-0.01 



Figure 5. The critical function ?/)2 calculated with the primes 



It is interesting to study the single contribution of a prime to the critical function 
^2- In Figure 6 we computed the contributions of the 10th prime {p = 29), of the 
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50th prime {p — 229) and of the 100th prime {p — 541), all the calculations were 
performed until q — 100. The computations indicate that not only the contributions 
decrease with increasing p but also that great primes have an influence only on big 
values of k. 




Figure 6. The contribution to the critical function ■02 of the 
primes p = 29 [black], p = 229 [red] and p = 541 [green] 



5. Infinite (3 limit 



In a numerical context we are also interested in the case of large P values. We 
start with the equation (7) in [16], given by: 
(5.1) 



(is) 



-Ijs ^p{s~p) ^ 2nf sH 



i{s + 2n) \ C(0) 



- 1 



and we set C — 



1 = ln27r — 1. Then, using the formula ^ = 
above (91(A) > 0!), we obtain: 



°° e'^'^dX as 



(5.2) 
where 

(5.3) 4 



fc=0 



1) 



1 



2n 



We now analyze ?A(fc), the absolute value of the critical function, at large (3 values 



where \V{^) ^ is valid. 



(5.4) V(fc) 



k « 



El , "-f 
— -k p 
„ P(a-p) 



E 



^ 2n{oL + 2n) 



g + 2n C , _g_ 
fc f _| /3 



Here the second and third term in the bracket converge for all ct > (in particular 
for I < CT < 1). If we choose a = 1, (5.4) would become in the /3 limit (supposing 
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that this limit may be performed and has a meaning): 
(5.5) 



lim ip{k) 

/3 — > oo 



(5.1) 

= hm 



i|lnC(^) + ^ 



= 7 



where 7 = 0.577216 is the Euler constant (see also [17]). 

If such a limit is permitted our conjecture is that for $H(s) > a + 5,5 > Q, as 
(3 ^00: 



(5.6) 



1/(^)1 



Bt 



-7 



where B is some constant and t = J(s). 

Since from the definition Pfc(^ + 1) = ^^=^iPfc_i(^) we obtain: 



fe=i 



la— s ,s — a, 
-— P,_i(— ) 



Then applying the Baez-Duarte inequality, i.e. |Pfc_i(z)| < 
for 9l(s) >a + 5,5>0 that: 



we have 



P 



i/(«)i~E 



00 B \ — 



fc=2 



+ 1 



V'(A;) ~ P 



a - s 



-V(fc) 



and finally: 
1/(^)1 



„|a-s|(A:-l) " ,,,, „|a-s|/3,,,. 



BJ(a-a)"+t2 



A similar (of course not rigorous) limit is formally obtained for ip2{k) using the 
primes along the lines for (3.4) to (3.10), which, as jS^oo is given by: 



<^ , 11 

lim k '3 dk= . 

00 a a — 1 a 



1 °° 1 

pprime g=l 



and thus [17]: 



lim 

a-> 1+ 



a 1 a - 1 ^ '^'P 2^ pocq I 

\ p prime g=l / 



lim - log((a - 1) C(a)) = 7 



We carried out some numerical experiments restricted to large fi values (until /3 = 
10^), using the first 3600 known zeros [18]. The computations in Figuro 7 indicato 
that for a fixed k, within the limit of accuracy of our computations, the difference 
between (5.4) and 7 approximately stabilizes to less than 0.001 indipendently from 
the choice of k. The difference is largely due only to the term involving the non- 
trivial zeros. That is if we need a higher precision we have to consider more non- 
trivial zeros in (5.4). 
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0.578 
0.5775 

0.577 
0.5765 




Figure 7. The convergence to 7 of ^p{k) for k = e^^ 
6. Conclusions 

In this work we have found some new representations of functions related to the 
Riemann Zeta function in terms of the Pochhammer's polynomials, i.e. for the Zeta 
function via the alternating series, for (1 — 2^~*)C(s), for ln((l — 2^~*)C(s)) and for 
the derivative of lii((s — l)C(s)). 

(1) A numerical experiment for the first function give satisfactory results both 
for the real part as well for the imaginary part even on the critical line 
9^(s) = i (we have used the values a = \, (3 = i and t up to 3{s) = t < 35). 

(2) In a formal limit of our representations (2.6) for the special case a = (3 = 2 
we obtain the Maslanka's representation of (s — l)C(s). 

(3) For the expansion of the derivative of the function ln((s — l)(^(s)) in terms 
of the Pochhammer's polynomials Pfe(s) we have found two expressions 
(■^1 and il)2) for the so called critical function: ipi in terms of the primes 
and il)2 in terms of the trivial as well as the non-trivial zeros. We have 
then carried out a numerical experiment which gives a very satisfactory 
agreements between the two, which up to very high values of k remain 
bounded. The existence of absolut upper bounds for the critical functions 
at A;-infinity may be considered as being equivalent to the truth of the RH. 

(4) Concerning the critical function in the large limit, using a = 1, we may 
conjecture that ^ time the derivative of ln((,s — l)C(s)), using the inequality 
for the Pochhammer's polynomials has, for fH(s) > 5 + (5, a bound of the 
form ^7 where 7 is the Euler constant and t = J(s). 
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